The modulated spin liquid: a new paradigm for URu 2 Si2 
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We argue that near a Kondo breakdown critical point, a spin liquid with spatial modulations can 
form. Unlike its uniform counterpart, we find that this occurs via a second order phase transition. 
The amount of entropy quenched when ordering is of the same magnitude as for an antiferromagnet. 
Moreover, the two states are competitive, and at low temperatures are separated by a first order 
phase transition. The modulated spin liquid we find breaks Z4 symmetry, as recently seen in the 
hidden order phase of URu2Si2. Based on this, we suggest that the modulated spin liquid is a viable 
candidate for this unique phase of matter. 



The hidden order (HO) phase in URu 2 Si2 is a long 
standing problem in condensed matter physics [T], and 
has recently received renewed attention from both ex- 
periment and theory. This metal can be characterized as 
a moderate heavy fermion with a Sommerfeld coefficient 
of ~ 180 mJ/mol K 2 . It undergoes a mean field- like sec- 
ond order phase transition at Tq = 17K, characterized 
by a large jump in the specific heat. The amount of en- 
tropy quenched at the transition is substantial (w 1.38 
J/mol K), which is of order of 24% of the entropy of a 
local / doublet. Despite many years of intense investi- 
gation, the nature of the hidden order remains contro- 
versial. A number of theoretical propositions have been 
made, which can be divided into itinerant types, where 
the order parameter originates from delocalized / elec- 
trons, and localized types, where it is believed that the 
local levels of a U ion are involved [H [3] . 

Experimental data show evidence for both itinerant 
and localized character of the order parameter. Under 
pressure, the phase diagram evolves from the hidden or- 
der phase to an antiferromagnet (AF) [4] . Recent inelas- 
tic neutron scattering (INS) measurements show that in 
the HO phase, a resonance occurs at a commensurate 
wave vector Qo = (1,0,0) [SJ|S], which transforms into 
a strong elastic AF signal for pressures P > 5 kbar. An 
inelastic resonance at Q* = (1 ± 0.4,0,0) occurs in both 
the HO and AF phases. The compensated nature of this 
metal leads to several quasi-nested portions of the Fermi 
surface, which could account for the Q vectors of these 
resonances, as well as the formation of an SDW [7J [5] . 
Recent STM experiments [9l [10] reveal the opening of 
a gap in the dl/dV characteristic at w To inside the 
HO phase. Moreover, careful analyses of specific heat [5] 
as well as thermal transport [TT] have revealed a strong 
analogy with the superconductor CeCoIns. This body 
of observations, combined with the opening of a Fermi 
surface gap inferred from conductivity [321 H3] , angle re- 
solved photoemission [14] , infrared spectroscopy [HI [16] , 
and the Hall effect [TTJ , constitute evidence for an itiner- 



ant mechanism. The localized viewpoint is based on the 
observation that under pressure, the first order transition 
between the HO and AF phases shows no distinct changes 
in the transport properties. The resistivity, for example, 
is continuous through the transition |18j . From this ob- 
servation, one can expect the HO phase to have strong 
similarities with the AF phase, in particular a doubling 
of the unit cell. 

In this paper, we offer a new idea for explaining the 
mysterious hidden order, which naturally interpolates be- 
tween the itinerant and localized viewpoints. The idea 
amounts to the observation that a spin liquid with spa- 
tial modulations that break translational symmetry can 
form via a second order transition with a large jump in 
the specific heat, while remaining "hidden" to most ex- 
perimental probes. 

The concept of a spin liquid dates from the early work 
of Fazekas and Anderson that a resonating valence bond 
(RVB) state might describe frustrated spin systems on a 
triangular lattice [TS] . Subsequent work by Anderson ex- 
tended this concept to describe high temperature cuprate 
superconductors [20] , In a seminal paper |21j . the uni- 
form spin liquid was understood as being a superposi- 
tion of bonding and anti-bonding valence bond singlets, 
stabilized by quantum fluctuations. This concept has 
re-emerged recently in the context of quantum critical 
points (QCP) in heavy fermions [22TI24] . The common 
understanding is that the anomalous properties of these 
metals arise from a competition between the formation 
of magnetic singlets between localized spins, and the for- 
mation of Kondo singlets. Recently, however, a body 
of experimental observations combined with theoretical 
insights has concluded that the formation of magnetic 
singlets is most probably the dominant mechanism for 
quenching the entropy of the local / spins [25] . The com- 
plexity of actinide ions, under the combined influence of 
spin-orbit, Hunds rules, and geometric frustration, cre- 
ate the optimal grounds to favor the formation of valence 
bond singlets that eventually quench the high tempera- 
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ture entropy . A spin liquid can simply be re-defined as a 
regime of localized spins in which the entropy is quenched 
despite the absence of long range order. 

In the RVB formalism, the spin liquid can be described 
within the t — J model [2T], which in the limit of half 
filling reduces to the Heisenberg model, with fermionic 
spin operators 
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where \\ a (xia) are creation (annihilation) operators for 
fermions with spin 1/2, charge zero, and gauge charge +e, 
commonly dubbed 'spinons'. For simplicity, we restrict 
ourselves to the square lattice in two dimensions, leaving 
a three dimensional generalization more appropriate for 
URu2Si2 to later work. Here er = ± is the spin index of 
the SU(2) representation, and spinons are subject to the 
constraint of one per site: xlaXia = L In the mean 
field approximation, the interaction term in ([!]) can be 
decoupled in a variety of effective fields to minimize the 
free energy, including valence bond singlets, Ylo-ixlo-Xjcr), 
valence bond "pairs", (xlaXja - Xig-Xja), or AF order, 
S ct /3(xL cr a/9X^e 4Q ' r '), with er the SU(2) spin matrix, 
and Q the AF wave vector. 

In the RVB theory for cuprate superconductors, 
uniform valence bond order parameters (XiaXja) — 
fo^ri^j+z where z is the index of nearest or next nearest 



neighbors, or flux phases (xlaXjc 
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*/ 4 ™ij where 



$ is the flux per plaquette and fiy = (r.; — rj)/a (with 
a the lattice spacing) is a bond orienting number, are 
commonly introduced. In the 7r-flux phase, for example, 
when a spinon cycles around a plaquette, the circulation 
of the phase of the order parameter generates a magnetic 
flux of magnitude n (Fig. [I]). Here we introduce another 
kind of valence bond order parameter, with real space 
modulations of the bond centers, which we denote as the 
modulated spin liquid (MSL): 
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We see that the value of the bond acquires an oscillat- 
ing sign from site to site. Since the phase on the bond 
is not oriented, no flux is generated when a spinon cy- 
cles around a plaquette. This order parameter doubles 
the unit cell associated with the dual lattice of bonds, 
breaking the Z4 symmetry of the underlying square lat- 
tice (Fig. [l]) . Whereas the flux phase is typical of lattices 
in two dimensions, the MSL can easily be generalized to 
three dimensions. 

The interaction term in ([I]) can be decoupled as either 
a spin liquid (SL) or an AF. For given J = Jsl + Jaf we 
consider J5 j, and Jaf as tuning parameters in our model, 
which thus acquire a phenomenological character. Con- 
sidering that the spin liquid can have both a uniform and 



a modulated component, the corresponding decouplings 
are performed via Hubbard-Stratonovich transformations 
on each bond ij, leading to the following Lagrangian 
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Here ifij is the Hubbard-Stratonovich field introduced 
for SL decoupling of the bond ij, and arises from 
the AF decoupling of the site i. In the following, these 
will be replaced by their constant, self-consistent, mean- 
field expressions, ipij = ~JsLj2a(xlaXja), and mi = 
JAFj2a a (XiaXicr)- Note that the SL field is defined on 
the dual lattice, whilst the AF one is defined on the ini- 
tial (square) lattice. In general, (fij can have a non-zero 
imaginary part, which is the case, for example, in a tt- 
flux phase. Here, we consider real SL fields only, which 
reflects a symmetry (pij = ipji on each bond. Such a sym- 
metry being incompatible with the occurrence of a mag- 
netic flux, the resulting modulated SL phase is necessarily 
of a different nature. Note that in two dimensions, flux 
phases might still coexist with the present MSL phase, 
but the study of this phenomenon is beyond the purpose 
of the present work. We introduce the Fourier trans- 
formed fields, ipq and m q . Our analysis concentrates onto 
the three following mean-field parameters: the uniform 
SL, 4>q = y(o,o)i the modulated SL, 0q = <^( W)7r ), and the 
Neel AF, Sq = m^^), noting that (w, it) has a lower free 
energy than (jr, 0). 

An intuitive description of the MSL phase, compared 
to other possible SL phases, can be obtained from 
finite size versions of our model. For two sites and 
Jaf = 0, the effective mean-field Hamiltonian ^ can be 
decomposed into bonding and anti-bonding eigenmodes, 

H12 = + tpi2 E CT (x A *XA* - XBaXBa) , with 

XBa = ^(Xla+X2a) and \Aa = ^ (Xla ~ Xia ) ■ 

These modes are reminiscent of the singlet, and 
the S z =0 triplet, which diagonalize the initial, two 
sites, Hamiltonian |l]). Invoking the saddle-point 
self-consistent relation for ip 12, we find a ground state 
energy, (i/12) = — j* L ■ The ±(^12 degeneracy reflects 
the U(l) local gauge symmetry of the Hamiltonian ([I]), 
which is invariant with respect to the transformation 
xla- ~x\a- The groundstate of the mean-field effective 
model can therefore arbitrarily be chosen as a bonding 
or an anti-bonding mode. Due to this gauge symmetry, 
for the two sites model, all SL mean-fields are equiva- 
lent and mimic the energy splitting between a singlet 
ground state and one of the triplet excited states. This 
equivalence does not hold anymore for bigger systems, 
where the number of sites (i.e., the number of local 
gauge symmetries) becomes smaller than the number of 
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(a) MSL (b) Flux Phase 
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FIG. 1: (Color online) the MSL phase (a) and the 7r-flux 
phase (b). Note the orientation of the bonds in the latter 
case. (c,d): the spinon dispersions for the MSL (bands 1 and 
2). Note the breaking of Z4 symmetry, and the small hole 
pockets for band 1 and the small electron pockets for band 2 
(white curves). The parameters are t'=0.1 and Jsl~5, with 
0o=2.11 and <£ Q =1.27. 

bonds. In the MSL case on the square lattice, however, 
the bonding and anti-bonding character oscillates from 
site to site, since the sign of the hopping parameter, tfij, 
is oscillating. 

The MSL can be considered as the true RVB parent of 
the AF phase. It can also be viewed as a liquid phase 
of dimers, where a resonance moves between different 
dimer coverings of the lattice. In three dimensions, the 
liquid phase of dimers can coexist with the breaking of 
translational invariance. 

Fig. [l] depicts the spinon dispersion of the MSL. The 
breaking of Z4 symmetry is obvious. Comparison of the 
free energy and specific heat jump at the transition re- 
veals that the MSL and the AF decouplings are almost 
degenerate at the mean field level. The MSL can thus be 
considered as the RVB parent of the AF order. 

The issue of whether a gauge invariant Lagrangian de- 
velops crossovers or phase transitions between the Higgs, 
confined and Coulomb phases is an old one and much 
related to the presence of instantons in the system [2"rj] . 
In the case of the MSL, we argue that the breaking of Z4 
symmetry is enough to ensure a second order transition 
|27j . This is to be contrasted with the 7r-flux phase or 
the possible condensation of the holons with modulations 
6 i e J< ^' ri [23]. In both of those cases, the order parameter 
would be sensitive to the effect of gauge fluctuations, and 
the presence of instantons at finite temperature generates 
a crossover. 

In the remainder of the paper, we discuss the compati- 
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FIG. 2: (Color online) (a) mean field phase diagram in 
the (T, Jaf) plane from Q for a square lattice model with 
t' — 0.1 and Jsl~5. Below the green curve (circles), <j>o turns 
on and below the blue one (squares), 4>q condenses. Sq con- 
denses below the red curve (diamonds). For low T this is a 
first order transition, at higher T a second order one. (b) 
Variation of <j>o, <pQ, and Sq as a function of Jaf for T=0.1. 

bility of the MSL with the HO phase in URU2S12, and its 
relation to the AF phase. To proceed, we solve a simpler 
model where the conduction electrons are ignored, and 
the f electrons are treated as a single orbital on a square 
lattice with Q — (tt, it). The resulting free energy per 
site can be written as: 

F = _^ EM1+e -^) + iM±fH + |L (4) 

where /3 = 1/ (/c^T), N is the number of sites, i runs over 
the two bands (each of which are Kramers degenerate) , 
and the spinon dispersions u>i(k) are given by 

Wi(fc) = y (e k + e k+Q ) - M/ (5) 
±\/f {£k £k+Q)2 + ^< N Q/^ 2 + WS Q 

with 

efe = —2(cos(k x a) +cos(k y a)) —At' cos(k x a) cos(k y a) (6) 

and /if the chemical potential of the / electrons adjusted 
so that the system is at half filling. Here, t' represents 
the next near neighbor contribution to <fo- For <j>Q, only 
near neighbor bonds are considered (with e"^ the first 
term of the previous equation), and to obtain a real field, 
the quantity 8 in ^ is set to tt/2. We minimize the free 
energy using Powell's method, and use a root finder to 
determine Hf for each choice of ((/>o,(I>q,Sq) [25] , 

In Fig. [T] we show the spinon dispersions for the MSL 
case. Note the pronounced breaking of Z4 symmetry, 
which has been recently detected in the HO phase from 
susceptibility measurements |29j . A large part of the 
Fermi surface (as defined when <pg is zero) is gapped 
upon ordering |T|, which gives a natural explanation for 
the amount of entropy quenched, and is consistent with 
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the Hall [T7] , thermal conductivity [TT] and quantum os- 
cillation [30] data that suggest that 90% of the carriers 
disappear at the transition. The dispersion for the AF 
phase is very similar, except in our simplistic approxi- 
mation, a full energy gap occurs once Sq > t'<j>o/2. Wc 
also find that the the order parameters in each phase 
have similar magnitude. One can tune between the two 
phases by varying Jaf relative to Jsl, as shown in Fig. [5] 
Note that (j) and 4>q are quasi-degenerate and condense 
at almost the same T. They would be equal if t' = 0. 
Above their condensation temperature, the free energy 
goes as — 2/csTln(2), which is just the free energy for 
a local / doublet. Sq condenses as Jaf increases, via 
a first order transition for low T, changing to a second 
order transition at higher T. In the AF phase, all three 
order parameters are at first non-zero, with 4>q eventually 
disappearing (at low T, we find a finite <fio). The compe- 
tition between the modulated component (4>q) and the 
AF order (Sq) that gives rise to the first order behavior 
is obvious from Fig. [2] This reproduces the qualitative 
features of the experimental phase diagram under pres- 
sure where a first order transition occurs between the HO 
and AF phases. Note that since both orders double the 
unit cell, no noticeable feature is expected to be seen in 
the electrical conductivity at this transition. 

From the more complete theory, we expect that the 
MSL phase will be stabilized near a Kondo breakdown 
QCP pT, reflecting the localization of the two 5f elec- 
trons per U site due to strong Coulomb forces. At the 
localization transition, the two 5f electrons decouple from 
the conduction electrons to form a spin liquid. The ef- 
fective hybridization is renormalized to zero at the QCP. 
The tunneling experiments confirm the opening of a hy- 
bridization gap at w To between heavy and light hole-like 
bands [9J [10] . We interpret the hybridization gap as the 
effective hybridization between the spinons and the con- 
duction electrons. Note that the spinons are very difficult 
to observe; they can be detected by transport only when 
they hybridize with the conduction electrons. This could 
explain why the heavy band seems to disappear above To 

INS experiments show an inelastic resonance at a com- 
mensurate wave vector Q = (1, 0, 0) in the HO phase [5J, 
which becomes an elastic peak in the AF phase. Our the- 
ory also has an inelastic resonance at Qaf since the AF 
state has a higher free energy in the HO phase. We spec- 
ulate that the incommensurate peaks are due to nesting 
of the spinon Fermi surface. 

In conclusion, the modulated spin liquid is an inter- 
esting phase of matter which has many properties com- 
patible with the hidden order phase in URu2Si2. The 
hidden nature of our order parameter simply relies on 
the fact that a spin liquid is hardly detectable if no long 
range symmetry is broken. We believe that it is a viable 
candidate for the solution of this long standing mystery. 
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